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ABSTRACT 



Aims. We investigate the size-density relation in extragalactic H n regions, with the aim of understanding the role of dust and different 
physical conditions in the ionized medium. 

Methods. First, we compiled several observational data sets for Galactic and extragalactic H n regions and confirm that extragalactic 
Hn regions follow the same size (f)-density (n) relation as Galactic ones (n oc D" 1 ), rather than a relation with constant luminosity 
(n oc D l 5 ). Motivated by the inability of static models to explain this, we then modelled the evolution of the size-density relation of 
H ii regions by considering their star formation history, the effects of dust, and pressure-driven expansion. The results are compared 
with our sample data whose size and density span roughly six orders of magnitude. 

Results. The extragalactic samples cannot be understood as an evolutionary sequence with a single initial condition. Thus, the size- 
density relation does not result from an evolutionary sequence of Hn regions but rather reflects a sequence with different initial gas 
densities ("density hierarchy"). We also find that the size of many Hn regions is limited by dust absorption of ionizing photons, rather 
than consumption by ionizing neutral hydrogen. Dust extinction of ionizing photons is particularly severe over the entire lifetime of 
compact H n regions with typical gas densities of ^ 10 3 cirr 3 . Hence, as long as the number of ionizing photons is used to trace massive 
star formation, much star-formation activity could be missed. Such compact dense environments, the ones most profoundly obscured 
by dust, have properties similar to "maximum-intensity starbursts". This implies that submillimeter and infrared wavelengths may be 
necessary to accurately assess star formation in these extreme conditions both locally and at high redshift. 

Key words, dust, extinction - galaxies: dwarf - galaxies: evolution - galaxies: ISM - galaxies: star clusters - H n regions 



1. Introduction 

Massive stars, young star clusters, and their associated H n re- 
gions are vital probes of recent star formation in nearby galax- 
ies and the distant universe. While Galactic Hn regions such as 
the Orion Nebula and RCW 49 are well studied because of their 
proximity, most Galactic work is still plagued by distance uncer- 
tainties (e.g., Anderson & Bania 2009), although the situation 
is improving rapidly (see, e.g., Hachisuka et al. 2006 Foster & 
MacWilliams 2006, Russeil et al. I2007I) . Hn regions in galax- 
ies in the Local Group can be studied in almost as much detail 
as those in the Galaxy, and the distance estimation is much less 
uncertain. However, even in Local Group galaxies, Super Star 
Clusters (SSCs), the most extreme examples of massive star for- 
mation, are absent. SSCs, with ~ 1O 5 -1O 6 M and Lyman con- 
tinuum photon rates N[ on ~ 10 52 — 10 53 s _1 enclosed in regions 
of $10 pc in radius, are generally not seen in quiescent environ- 
ments such as the Milky Way, M 31, and M 33. Even R136 in 30 
Doradus, the most massive Local Group star cluster (1O 4 5 M 
and 10 51 4 s _1 : Massey & Hunter 1998 ), falls short of the typical 
properties of SSCs. Hence, to study the wide range of manifes- 
tations of massive star formation, it is necessary to examine star 
clusters and Hn regions in galaxies beyond the Local Group. 
Of necessity, the price to be paid is detail; the advantage to be 
gained is the wide variety of Star-Forming (SF) complexes that 
can be studied. 

The total mass of a cluster of massive stars can be observa- 
tionally quantified by estimating the number of ionizing photons. 
A simple Stromgren-sphere argument (Eq.|2| would suggest that 



7Y ion oc D 3 n 2 , where Ni on is the number of ionizing photons emit- 
ted per unit time, D is the diameter of the H n region, and n e is 
the electron number density in the H n region. Thus, in princi- 
ple, by examining n e and D of H n regions, we can estimate the 
ionizing photon luminosity. However, observationally, the situ- 
ation is not straightforward. The diameter (D) and the electron 
density (n e ) of Galactic H n regions are known to be negatively 
correlated with a roughly unit slope: « e oc D l (Garay & Lizano 
[19991 Ki m & Koo |2"00Tl Martin-Hernandez et al. |2003t Dopita et 
al. l2006b . This would not be expected from a simple Stromgren- 
sphere argument with A^i on constant: n e cc D 15 . Various possible 
explanations for the observed shallower slope are presented in 
Martin-Hernandez et al. (2003 ), including optical depth effects, 
dumpiness, dust extinction, and stellar content. Dust mixed in 
with the ionized gas has also been proposed as an explanation 
by Arthur et al. (2004, hereafter A04) who argue that significant 
absorption of ionizing photons by dust grains in the densest H n 
regions flattens the size-density relation. 

Previous work on extragalactic Hn regions suggests that 
they follow the same size-density correlation as Galactic ones. 
Kennicutt (1984) found that Hn regions in spiral disks follow 
a n e oc D relation, but extend the Galactic trend to lower 
densities and larger sizes. More recently, Gilbert & Graham 
(2007) examined the SSCs in the Antennae galaxies, a proto- 
typical starburst merger, and found that these also follow the 
« e oc D l relation. However, the SSCs in the Antennae have den- 
sities (~ 40-400crrT 3 ) and sizes (25-100 pc) that place them on 
a different location in the n e -D plane than the Kennicutt sample. 
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In fact, Gilbert & Graham ( 120071 1 conclude that the SSCs in the 
Antennae constitute a new class of massive H n regions that is 
distinct from the Galactic and typical extragalactic population. 

In this paper, we examine the size-density relation in ex- 
tragalactic H ii regions, and explore several mechanisms which 
have been proposed to explain it, including dumpiness, dust 
extinction, and stellar content. In particular, we focus on the 
Hn regions in Blue Compact Dwarf galaxies (BCDs); follow- 
ing Gilbert & Graham (2007 ), we shall refer to such regions as 
Emission-Line Clusters (ELCs). In fact, most of the current star 
formation in many BCDs occurs in such ELCs, rather than in 
the underlying diffuse component. Previously we studied a more 
limited sample, and found that the H n regions in BCDs also fol- 
low the same size-density relation namely, n e oc D 1 (Hunt et al. 
2003 ). Here we triple the sample size used in our previous work, 
and compare the BCD ELCs with other types of extragalactic 
H ii regions, such as those in spiral disks and known SSCs. 

Focusing on BCDs enables us to study another aspect of the 
formation of star clusters. Because they are generally metal poor, 
with oxygen abundances [12 + log(0/H)] ranging from 7.2 to 
8.5 (Izotov, Thuan, & Stasiriska |20U71 ), BCDs provide a link be- 
tween high-redshift metal-free primeval galaxies and the metal- 
enriched SF galaxies in the local universe. Because metallicity 
is not expected to be the only driver of star-formation and clus- 
ter properties, we can compare different metallicities together 
with other parameters, and better disentangle the effects of metal 
abundance. 

In fact, the process of ionization itself in metal-poor objects 
is also important in the cosmological context. Because they are 
relatively chemically unevolved, Hn regions in BCDs could en- 
able us to infer some characteristics of the ionization processes at 
high z. Since the typical virial temperature of the first-generation 
objects in the Universe is % 10 4 K (e.g., Tegmark et al. 1997; 
Yoshida et al. 2003 ), the photoionization which raises the gas 
temperature to ~ 10 4 K prohibits the gas from collapsing to form 
stars (Omukai & Nishi 1999 ). The gas density structure is mod- 
ified by the pressure-driven expansion of H n regions (Kitayama 
et al. 120041 ). Such effects on the gas structure are of fundamen- 
tal importance for considering the subsequent star formation and 
the reionization of the Universe. 

The relation n e oc D 1 for Galactic and extragalactic Hn re- 
gions can be interpreted as a constant ionized-gas column den- 
sity. However, the physical mechanism that produces the con- 
stant column density is not clear, nor is it clear why ELCs in ex- 
ternal galaxies and Galactic H n regions should follow a similar 
relation with column density. It is also not understood whether 
or not there are different classes of H ii regions which would oc- 
cupy distinct zones of parameter space in the size-density plane, 
and if there are, how they might be related. 

The aim of this paper is to understand the size-density re- 
lation of extragalactic Hn regions, and place it in the context 
of star-cluster and ELC formation. The paper is organized as 
follows. First, in Sect. |2j we describe the observational sam- 
ples of extragalactic H n regions, together with the compilation 
of Galactic H n regions for comparison. In Sect. [3j we examine 
gross trends in the data, which we will later interpret in the light 
of our evolutionary models. We predict the size-density relation 
for static dusty H n regions by assuming a constant luminosity 
of the central sources in Sect. [4] Then, in Sect. [5] we extend 
the model to include pressure-driven expansion and star forma- 
tion history. Some basic results of this evolutionary model are 
given in Sect. [6] and compared with the extragalactic observa- 
tional data. In Sect. [7] we discuss the models and their impli- 



cations in various contexts. Finally we give our conclusions in 
Sect. [8] We adopt a Hubble constant of H = 73 km s _1 Mpc" 1 . 



2. The data 

We have assembled several samples of H n regions from the lit- 
erature, including two extragalactic and two Galactic radio data 
sets, and one extragalactic optical sample. Another extragalac- 
tic optical sample is presented here for the first time, with sizes 
measured from HST images and electron densities derived from 
emission measures calculated with published long-slit optical 
spectra. 

2.1. Extragalactic radio data sets 

One extragalactic radio data set comprises those galaxies with 
known "ultra-dense Hn regions" or "radio super nebulae" (see 
Kobulnicky & Johnson [19991 Beck et al. 120021 ). Many of these 
are found in low-luminosity low-metallicity BCDs, but some 
reside in metal-rich starbursts such as NGC253 and M82 and 
normal spiral disks, including NGC4214 and NGC6946. The 
common feature of these multifrequency radio continuum ob- 
servations is a rising spectrum at low frequencies, and a rela- 
tively flat spectrum at higher ones. This implies that the pre- 
dominant emission mechanism is thermal bremsstrahlung from 
ionized gas, and that the electron density and size are such that 
there is a turnover in the radio spectrum at some frequency v t . 
This "turnover frequency" v, corresponds to an optical depth r v 
of unity, and defines the frequency where the spectrum changes 
from optically thick to optically thin, v, depends on the Emission 
Measure (EM = J n 2 e d() and the electron temperature T e . The 
spectra of "classical" H n regions such as Orion turn over at very 
low frequencies, ~0.3GHz, while the compact dense regions in 
these galaxies have v, >5 GHz. All radio observations described 
below measure the high-frequency optically thin part of the radio 
spectrum (F Y oc v~ 01 ), as well as the transition region (v < t v ) 
toward lower frequencies where the spectrum becomes optically 
thick. 

In some cases, with high-resolution observations (IIZw40: 
Beck et al. 120021 NGC 5253: Turner et al. 120001 ), the radio emis- 
sion is resolved at 15 GHz. Thus, the authors were able to es- 
timate the size D of the emitting region by fitting a Gaussian, 
and from the observed flux, determine the EM, and infer the 
rms electron density <n e >. In other cases (NGC 4214, NGC 1741, 
Mr k8, M rk33, VIIZwl9, Pox 4, Tol35, Mrkl236: B eck et 
al. 120001 NGC 6946, NG C 253 , M33: Johns on et al. I200T1 
SBS 0335-052: Hunt et al. EjQl Johnson et al. 120091 , the spa- 
tial resolution is insufficient to resolve the regions. Then the ra- 
dio spectrum can be fit by models of homogeneous, isothermal, 
dust-free, ionization bounded regions of ionized gas, to obtain 
the turnover frequency v,, the EM, (n e ), and infer the size D of 
the region (e.g., Deeg et al. 1993; Johnson et al. 2001; Hunt et 
al. 120041 ). Alternatively, the optically thick and optically thin re- 
gions of the radio spectrum can be separated to constrain v,, and 
thus infer the emission measure EM, size D, and rms electron 
density <n e > (e.g., Gordon [19881 Beck et al. 120001 ). When the 
size is not directly measured, that is to say when the sources 
are unresolved, the authors estimate sizes and densities from 
fits of multi-frequency radio continuum spectra; these are con- 
sequently not independent parameters, but rather negatively cor- 
related. Nevertheless, the logarithmic slope between density and 
size expected from this degeneracy would be - 1 .5 which is sig- 
nificantly steeper than that observed (see below). The sizes D 
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and densities (n e ) for NGC5253 and He 2-10 have been mod- 
eled also from radio recombination line observations, and are 
consistent with those inferred from continuum fitting (Mohan et 
al. l200Tb . 

There is one galaxy in our data set, IZw 18, in which multi- 
frequency observations show no sign of a rising spectrum (Hunt 
et al. 120051 1. Here also the data have been fit to a model of an 
homogeneous, isothermal, dust-free ionization bounded region 
of ionized gas, as described above. 

We will refer to this sample of (except for IZw 18) rising- 
spectrum sources as the "radio sample" (16 galaxies); its mean 
oxygen abundance is 12+log(0/H) = 8.24+0.5, or -0.22 Zefl 
When there are multiple observations for a single object, we 
usually list these data as different data points. However, M 33, 
NGC253, and NGC6946 contain several candidates for ultra- 
compact H ii regions, but we show each galaxy as a single aver- 
age according to the figures in Johnson et al. (2001). The data 
for the radio sample are reported in Table [T] all sizes have been 
corrected to the distance scale used here. 

The other extragalactic radio sample comprises H n regions 
in the Small and Large Magellanic Clouds (SMC and LMC, 
respectively) (Martin-Hernandez et al. 2005), and in the super- 
giant Hn region NGC604 in M33 (Church well & GossQll9]>- 
As with the previous data set, these are also radio continuum 
observations, but at a single frequency, 5 GHz (SMC/LMC) or 
8.4 GHz (NGC 604/M 33). We adopt the results given in the orig- 
inal papers for sizes D and densities (n e ). The regions in these 
Local Group galaxies are resolved, and D is measured from 
high-resolution interferometric maps by fitting two-dimensional 
Gaussians, including beam deconvolution. The densities are in- 
ferred from the observed flux densities by assuming that all radio 
emission is optically thin bremsstrahlung, arising in a dust-free, 
ionization bounded, homogeneous region with size D (Mezger 
& Henderson 1967J. This sample will be called the "Local- 
Group sample"; with M33, the LMC, and the SMC, its mean 
oxygen abundance is 12+log(0/H) = 8.41 ±0.3, corresponding to 
0.33 Zq. 

2.2. Galactic radio samples 

We have also included two size-density data sets of Galactic 
Hn regions as comparison samples. The first Galactic sample is 
taken from the compilation by Garay & Lizano ( 1999), and con- 
sists of compact H n regions observed with high angular resolu- 
tion in either the H66a or H76a lines, as well as in the radio con- 
tinuum. The second sample is a set of (ultra) compact Galactic 
Hn regions, observed in the 21 cm radio continuum by Kim & 
Koo (120011 1. For both samples, we have adopted the authors' H n 
region parameters; the electron densities were derived assum- 
ing optically thin, dust-free, homogeneous, ionization-bounded 
nebulae (e.g., Mezger & Henderson 1967), and the sizes were 
directly measured from interferometric radio images by fitting 
Gaussians, deconvolved with the beam size. As in the Local- 
Group sample, such measurements obviate the potential density- 
size degeneracy with logarithmic slope -1.5 that could arise 
were the sizes not measured independently. 

2.3. HST extragalactic optical sample 

The main optical sample includes those star-forming dwarf 
galaxies with optical spectra and usable high-resolution Hubble 

1 In this and future relative solar abundances, we will adopt the solar 
calibration of Anders & Grevesse dl989t . 



Space Telescope (HST) archival data, obtained either with the 
Advanced Camera for Surveys (ACS), the Wide Field Planetary 
Camera 2 (WFPC2), or in the near-infrared with the NICMOS 
array. Most of the images were retrieved from the Hubble Legacy 
Archival We obtained images for 26 galaxies, but use only 
23 of them; 3 had only F160W images, and when images at 
other wavelengths were available for comparison, the F160W 
images gave consistently larger sizes than the other wavelengths. 
Virtually all of the galaxies are BCDs. 

The high resolution of HST is crucial in order to resolve the 
ELCs in the sample objects. Most of the galaxies are dominated 
by a single bright SF complex; because we want to compare with 
ground-based long-slit optical spectra, this is the region we fo- 
cus on, rather than examining the more diffuse emission. There 
were generally no Ha images available, so we were forced to 
use the continuum to determine the size of the SF region. Since 
we wanted to match the size measurement as much as possi- 
ble with the spectroscopic slit, we adopted a one-dimensional 
method rather than two-dimensional models as in the radio. For 
each galaxy, we measured the linear extent by fitting the surface 
brightness profiles in two orthogonal cuts with Lorentzian and 
Gaussian profiles. Although Lorentzians fit the extended wings 
of the profiles better than Gaussians, they give widths that are 
-8% smaller. Hence, we adopted the Gaussian fits to be com- 
patible with the analogous fits for the radio source sizes. The 
HST profiles are extracted in 1-arcsec wide rectangular aper- 
tures, which match the spectroscopic slit (i.e., with roughly the 
same spatial resolution as the H/3 measurements adopted below). 
The diameter of the region is defined as the geometrical mean 
of the full widths at half maximum (FWHMs) of these two or- 
thogonal (Gaussian) profiles. If data are available in more than 
one band, we adopt the longest-wavelength data with the highest 
spatial resolution, although there is no significant trend of sizes 
with filter band (except for F160W as noted above). Despite 
the high spatial resolution of the HST, we were unable to re- 
solve the brightest complex in a few distant galaxies with par- 
ticularly compact SF regions. In those cases, we will be overes- 
timating the region size, and consequently underestimating the 
root-mean-square electron density; these will be discarded in the 
analysis. 

In general, it should be emphasized that the size measure- 
ment is a delicate and difficult procedure. Many of the objects 
have blended regions which HST did not resolve, but which 
would fall within a ground-based spectroscopic slit. Moreover, 
because most of the spectroscopic observations did not give the 
position angle of the slit, we had to use subjective judgment to 
determine the angles for our virtual cut apertures. It is also true 
that ionized gas in an ELC tends to be more extended than the 
underlying nebular continuum and stellar emission, so we are 
probably underestimating the size with our method and account 
for this empirically (see below). All these considerations make 
the diameter determinations only good to a factor of 2 or so, but 
the consistency with the other samples lends confidence to the 
procedure. 

The root-mean-square (rms) number densities (n e ) are cal- 
culated from long-slit observations of the optical H/3 recombina- 
tion line. Following Kennicutt ( 1984), we convert the H/3 surface 
brightness over the spectroscopic slit area to volume emission 
measure. Hydrogen emissivities are calculated with ionized-gas 



2 The Hubble Legacy Archive is a collaboration between the 
Space Telescope Science Institute (STScI/NASA), the Space Telescope 
European Coordinating Facility (ST-ECF/ESA) and the Canadian 
Astronomy Data Centre (CADC/NRC/CSA). 
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temperatures inferred from optical emission lines, as given by 
published tables (see Table Extinction is corrected for with 
published values of c(HyS), and ionized helium with a multiplica- 
tive factor of 1.08. The continuum size of an ELC is generally 
1 .5 to 2 times smaller than the ionized gas extent (Tenorio-Tagle 
et al. 120061 Silich et al. 120071) . so to convert the volume emis- 
sion measure to rms density (re e ), we have adopted a region size 
1 .5 times as large as actually measured in the continuum. To ac- 
count for the larger extension of the ionized gas relative to the 
continuum, we also used this enlarged size as the "true size" of 
the Hn regions measured from the images. Because of the con- 
siderable uncertainties in this entire procedure, the densities (n e ) 
are probably only good to roughly a factor of 2, being slightly 
less uncertain than the diameters because of the square-root de- 
pendence on EM. Nevertheless, this is similar to the uncertainty 
in the rising-spectrum radio sample, and to the comparison op- 
tical sample described in Sect. 12.41 The sample, which we call 
the "HST sample", is given in Table [2] together with the refer- 
ences for the spectroscopic data. The mean oxygen abundance 
of the HST sample (23 galaxies) is 12+log(0/H) = 7.89±0.35, 
-0.09 Z . 

The rms densities (« e ) for the HST sample are plotted against 
the densities determined from the [S n] optical emission lines 
n e ([S ii]) in Fig. Q] Similarly to previous work (Kennicutt 1984; 
Rozas et al. 1 19981 1, the densities inferred from the [S n] line ratio 
n e ([S ii]) are much higher than the rms values (« e ) inferred from 
the emission measure. This is because the densities measured 
in situ from line ratios tend to be weighted toward high-density 
high surface-brightness knots which occupy a small fraction of 
the total volume (e.g., Zaritsky et al. |1994| Kennicutt [1984| ). The 
difference between the two kinds of measurements suggests that 
the ionized gas is clumpy, with dense knots in a more diffuse en- 
velope (Kennicutt |1984t Zaritsky et al. |1994t Rozas et al. |1998D . 
In a homogeneous medium with optically thin dense clumps, the 
volume filling factor (FF) relates the rms density and the sulfur- 
derived one: yf5 = {n e )/n e ([S ii]), where 6 signifies the filling 
factor. Constant volume FFs are shown in Fig.[T] the data appear 
equally distributed from FFs of roughly unity to 10" 3 . The fill- 
ing factors are slightly lower, although comparable to those in 
Hn regions in quiescent spiral disks (Kennicutt 1984). In gen- 
eral, the HST sample follows the same trends in size and density 
as the other samples, and appears to be consistent with them. 
The six galaxies with FWHM <4.5 pixels in the HST images 
are marked with an arrow in Fig.Q] We are underestimating the 
rms densities and overestimating the sizes for these unresolved 
sources, and they are not considered in subsequent analysis. 



2.4. Comparison extragalactic optical sample 

The second optical sample is taken from the cornerstone 
study of giant H n regions in nearby spiral galaxies (Kennicutt 
1984). With ground-based photographic Ha emission-line im- 
ages, Kennicutt used a spherically symmetric shell model and 
solved the Abelian integral for the emission-measure profile. 
Most of the profiles are monotonically decreasing with ra- 
dius. One of the ELCs measured by Kennicutt, M82-A, has 
also been recently measured by another group; the old val- 
ues of (450 pc, 16 cm 3 ; Kennicutt 119841 ) are now found to be 
(4.5 pc, 1800 cm 3 ; Silich et al. 120071 ) . Another object is in com- 
mon with our sample, Mrk71 (NGC2366A=NGC2363): the 
old values are (560 pc, 4 cm" 3 ), in contrast with our new esti- 
mate of 14.4 pc (this includes the doubling described above), 
149 cm" 3 . The HST image of Mrk71 gives a diameter of 



1000 - 




10 100 1000 

n e ([SII]) [cm-*] 

Fig. 1. RMS densities (n e ) vs. published densities inferred from 
the [S n] lines at /U6717, 6731 A for the HST sample. The di- 
agonal lines correspond to constant volume filling factors of 1 
(dotted), 0.1 (short dashed), 0.01 (long dashed), and 0.001 (dot- 
dashed). The filled circles correspond to IZwl8 (n e ([Sn]) = 
120 cm" 3 ), and SBS 0335-052 (ra e ([Sll]) = 412 cm" 3 ), respec- 
tively. The red circles correspond to n e ([Sn]) >125cm" 3 , and 
the blue ones to n e ([Sn]) <125cm" 3 . The arrows show those 
galaxies with unresolved sizes, and thus underestimated (n e ); 
these objects are not considered further. 



7.5 pixels (~0'.'75), and clearly resolves the SF complex; this cor- 
responds to only the brightest portion of the region measured by 
Kennicutt, apparently not resolved by the ground-based photo- 
graphic images. It is clear that the larger the region over which 
the density is averaged, the smaller the rms density. It is notewor- 
thy that all these measurements, old and new, follow the same 
relation between size and density, as discussed below. 



3. Overall empirical trends 

The size-density relation of the samples is shown in Fig. [2] All 
Galactic and extragalactic samples can be fit by n e x D~ l to 
within the uncertainties in the slope. This behavior was already 
noted for the Galactic samples by Garay & Lizano ( 119991 ) and 
by Kim & Koo (1200 lj . Fig. [2] shows the best-fit regressions for 
the Kim & Koo sample, the Garay & Lizano sample, and the ex- 
tragalactic radio sample, but with the slopes slightly tweaked to 
be exactly unity; the intercepts for the three regressions (logn e 
[cm -3 ] at diameter = 1 pc) are 2.8, 3.5, and 4.4, from bottom to 
top, respectively. At face value, these offsets would imply that 
the mean column densities in the ionized gas increase by almost 
two orders of magnitude, going from compact to ultra-compact 
Galactic Hn regions, to ultra-dense extragalactic Hn regions 
(the radio sample). The size-density relation in all Hn regions, 
Galactic and extragalactic, is clearly flatter than that produced in 
a homogeneous Stromgren sphere with a constant luminosity of 
ionizing photons (n e oc D 15 ; Sect. [4]). 

The trends of the different samples also suggest that metal- 
licity is not the key factor in H n region properties. While at ap- 
proximately the same metal abundance, some of the regions in 
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Fig. 2. Densities n e vs. size (diameter) of the six Hn-region 
samples. The extragalactic radio sample (except for M 33 and 
IZw 18) is shown by filled (red) triangles; (the radio data for) 
IZw 18 by a filled blue triangle; M 33 by open (blue) squares; the 
Garay & Lizano d 19991 1 sample by ( green) x; the LMC/SMC by 
(black) open stars; the Kim & Koo (120011 ) sample by (cyan) +; 
the Kennicutt (1984) sample by open (cyan) triangles; the HST 
sample by open (red and blue) circles, except for IZw 18 which 
is shown by a filled (blue) circle. SBS 0035-052 has been ex- 
cluded from our HST sample because it is unresolved. The large 
rectangle shows the range of sizes and densities for the SSCs 
in the Antennae galaxies (Gilbert & Graham 2007). Root-mean- 
square densities (« e ) ( cm" 3 ) are shown for the HST sample, with 
sizes and (n e ) calculated as described in the text. The red cir- 
cles have [S n] -inferred densities of >125 cm 3 , and blue circles 
<125 crrr 3 . The regression lines are unit slopes (the same as the 
best fit, to within the errors) with intercepts of 2.8, 3.5, and 4.4. 
A factor-of-two error bar is shown in the lower left corner. 



the Local-Group sample lie closer to the dense compact Galactic 
sample, while others are coincident with the less dense Galactic 
sample. Both Galactic samples are of roughly solar metallicity, 
but their locations differ from one another in the size-density 
plane. Lastly, the mean metallicity of the ultra-dense radio sam- 
ple is only slightly lower than solar, but they lie far away from 
the locus of the (roughly solar abundance) H n regions in spiral 
disks. 

Except for the radio sample, most of the extragalactic H n re- 
gions follow the same size-density trend as the Galactic ones. In 
particular, a large part of the HST sample BCDs are located at the 
extension of the size-density relation of the Galactic and Local- 
Group samples, coincident with the H n regions found in spiral 
disks. The correlation of size and density in the HST sample 
alone is highly significant; with a parametric correlation coeffi- 
cient r - -0.66, the (one-tailed) significance level is >99.8%. 



3.1. Emission measure and density systematics 

All the densities in the size-density relation presented here are 
derived from the EM of either free-free radio emission or hy- 
drogen recombination lines in the optical. Because EMoc D, we 
might expect the size-density relation to result from constant 
EM over a sample, combined with a constant luminosity as in 
the Stromgren argument. The first would result in (« e ) °c D -1 ^ 2 , 
and the second would give (n e )oc £T 3 ^ 2 ; combining the effects 
would tend to flatten the slope from the Stromgren relation. 

However, we can exclude this as the reason for the unit 
slope in the data presented above, and in previous work by other 
groups. The EMs in the Galactic radio samples vary by more 
than four orders of magnitude; the same is true for the optically- 
inferred EMs in the HST sample. This would refute the hypoth- 
esis of a constant EM. Moreover, we have verified that the den- 
sities inferred from the S n lines for the HST sample are also 
correlated with the size D. These optical measurements are in- 
dependent of the EM inferred from the hydrogen lines, and thus 
should provide a robust check of systematics. Because the sul- 
fur lines do not probe densities significantly below ~ 100 cm" 3 , 
we exclude Sn densities with values <50cirT 3 , and find a cor- 
relation coefficient of r — -0.46. This is a weaker correlation 
than the one with rms densities (n e ), but still significant at the 
96% level. Since the S n densities are independent of the emis- 
sion measure, we conclude that the size-density relation is not 
spuriously induced by the method used to infer rms electron den- 
sities. 



3.2. Scale-free star formation 

The power-law size-density relation of H n regions suggests that 
massive star formation is self-similar, that is, there is no char- 
acteristic scale of star formation. This scale-free nature was al- 
ready noted by Larson (1981) who found an approximately unit 
slope between the rms H2 volume density and the size of molec- 
ular clouds: (re(H 2 )>(cm- 3 ) = 3400£(pc) -u . Kim & Koo d200TT ) 
found a similar relation relating Galactic H n region density and 
size, and argued that it reflects a variation of the ambient den- 
sity, rather than an evolutionary effect. Indeed, the similarity of 
the size-density trend for molecular clouds and H 11 regions sug- 
gests that the star-formation processes retain an imprint from the 
molecular environment in which they take place. 

This scale-free nature of H 11 regions is also supported by 
the data presented here. Three BCDs, He 2-10, SBS 0335-052, 
and IIZw40, host both ultra-dense radio nebulae (see Table [TJ, 
and optically-visible Hn regions (see Table |2). All these data 
follow the same size-density relation, but with different offsets. 
This implies that when the same regions are probed with longer 
dust-penetrating wavelengths and higher spatial resolution, they 
turn out to be denser and smaller, but with the same size-density 
relation as for the larger complexes. H 11 regions and hierarchical 
star formation will be discussed further in Sect. 17.41 



4. Static models 

First, we interpret the size-density relation of the Hn regions 
compiled in Sect. [2] by using simple theoretical arguments. In 
particular, we relate the size and density of H 11 regions for an 
ionizing point source embedded in a uniform and static medium 
with constant N lon . In such a situation, the radius of the ionized 
region can be estimated by the Stromgren radius (Spitzer 1978 ). 
We also include the effect of dust extinction, which is thought to 
be important in determining the size of H 11 regions (e.g., Inoue et 
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al. 120011 A04). The models described here are static models, and 
assume a number of ionizing photons Ni on constant over time. 

4. 1 . Size-density relation of dusty static H ii regions 

Here, we estimate the radius up to which the central source can 
ionize in a dusty uniform medium; this radius is called the ion- 
ization radius, r;. Before estimating r;, we define the Stromgren 
radius, rs, as 

4-7T 

—rln e n H a (2) = N ion , (1) 

where ff (2) = 2.6 x 10~ 13 cm 3 s is the recombination coeffi- 
cient excluding captures to the ground (n — 1) level (Case B), 
«h is the number density of hydrogen nuclei, and « e is the elec- 
tron density estimated by « e = 1.08«h, where the factor 1.08 
indicates the correction for helium ionization. Numerically, the 
Stromgren radius is estimated as 




In the absence of dust, the ionizing radius r; would be equal 
to rs. However, since a certain fraction of ionizing photons are 
absorbed by dust grains, the ionization radius r; is reduced by a 
factor of y;(< 1), which is determined by solving (Petrosian et 
al. H972l Spitzer fl978l A04) 

3 f 'y 2 e yTs "dy = 1, (3) 
Jo 

where rsd is the dust optical depth for the ionizing photons over 
a path length equal to the Stromgren radius. This integral can be 
expanded: 

1 = ^h e ™ - + 4~ (e Tsd - yi - 1) (4) 

rsd T Sd T Sd 

In the optically thin limit (rsa « 1), Eq. (0]) becomes yi - 1 - 
i"Sd/4. Thus, we confirm thaty; — > 1 (i.e., n — > rs) as Tsd — > 0. 

Once we obtain y; for a value of Tsd (given in Sect. 14.21 ), the 
ionization radius can be estimated as 

n = yi^s . (5) 

We also define the optical depth over the ionization radius as 
Tdi = yiTsd ■ (6) 

4.2. Dust optical depth rsd 

In the above, we have left rsd undetermined. Hirashita et al. 
(2001, see Appendix lAl for an alternative derivation) estimate 
it under a uniform dust-to-gas mass ratio D as 

The dust-to-gas ratio of the solar neighborhood is assumed to be 
2) = 6 10~ 3 (Spitzer 1978). We adopt various constant values 

3 This value is for a gas temperature of T = 10000 K (Spitzer 1978 ; 
Osterbrock & Ferland 2006). If we adopt T = 8000 K and 20000 K, 
we find a (2) = 3.1 x 10~ 13 and 1.4 x 10~ 13 cm 3 s" 1 , respectively. The 
temperature dependence of or (2) does not change our conclusions as long 
as we assume a typical gas temperature ~ 10000 K in Hn regions. 



for D and do not follow the time evolution of D in order to avoid 
uncertainty in the parameters concerning the chemical evolution 
models. 

From the equations in Sects. 14. 11 and |4~2l we can infer the 
qualitative behavior of r\ as a function of N mn . As Mon increases, 
rs and rsd increase (Eqs.[2]and[7]( with oc N i( L. Inspection of Eq. 
(01) suggests that if Tsd ^ 1, y\ drops in a very sensitive manner 
with an increase of Tsd because of the exponential dependence. 
Thus, if Tsd » 1, n increases only slightly with N- lon because 
yj decreases significantly. On the contrary, y\ ~ 1 if rsd < 1, 
since the dust extinction is not severe. In this case, r; is roughly 
proportional to Mem 1/3 ■ 

4.3. Dust-to-gas ratio, metallicity and filling factor 

It is important to realize that T> in Tsd (Eq. [7]) should not be 
interpreted strictly as a dust-to-gas ratio dependent on metal- 
licity. As implied by Fig. [1] the ionized gas must be clumpy, 
with concentrations of dense gas embedded in a more tenuous 
medium (Sect. [2731 . Following Kennicutt ( 1984 ) and Osterbrock 
& Flather (119591 ), we assume that the emission (and mass) of 
the ionized gas is dominated by the dense clumps. In this case, 
the much less dense inter-clump region would provide a neg- 
ligible contribution to the gas emission and mass. Implicit in 
this assumption is the optically-thin nature of the clumps (c.f., 
Giammanco et al. [2004). Therefore, the rms formulation of our 
models together with the assumption of a uniform distribution 
would dictate the introduction of a gas volume filling factor 6. 
Eq. (0 then becomes: 

— (M^w^T^T- «-> 

Because only a fraction of the volume 6 (6 1 ^ of the path length) 
is occupied by the dense gas, a factor of 6~ l t 6 must be intro- 
duced: r sd oc n^ nse (S 1/3 r s ) oc «™ s <r 1/2 <5 I/3 r s oc «™ s r s <T I/6 
("rms" and "dense" distinguish between the rms density and the 
density in the clumps). The rather mild dependence on 6 means 
that a filling factor of 0.001 would cause only an increase of a 
factor of 3 in rsd- 

For convenience, we define the dust-to-gas ratio + filling fac- 
tor normalized to the solar neighborhood value, k, as 

k = 6- 1/6 D/D Q . (9) 

In our models, we have assumed a uniform dust distribution 
within the H n region, since we have maintained a constant dust- 
to-gas ratio + filling factor. However, because of stellar winds or 
grain evaporation, the central cavity surrounding the star clus- 
ter could be devoid of dust (e.g., Natta & Panagia [19761 Inoue 
2002 ). In Galactic H n regions excited by single stars, volumes 
with sizes ~ 0.5 - 0.8 rs can be dust-free, which would imply 
dust filling factors as low as ~40%. Because of grain destruction 
effects and strong stellar winds, the dust-free volumes could be 
even larger in extragalactic H n regions ionized by massive com- 
pact star clusters with 10 3 - 10 4 O stars. In this case, we can still 
use our formulation, but should interpret the "effective" dust-to- 
gas ratio k as being small because of the central cavity. Thus, 
as mentioned above, the interpretation of small k is not unique: 
it could be due either to a low dust content, a relatively high 
gas filling factor, or to a non-uniform dust distribution caused 
perhaps by a dust cavity. It is likely that all three mechanisms 
contribute to the meaning of k in our formalism. 
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Fig. 3. Relation between the rms electron number density 
(n e ) and the ionization diameter D; for various samples. 
Symbols are as in Fig. [2] In addition to the observational 
data, some theoretical predictions for the static models are 
shown in each panel: (a) Left panel: the dotted (cyan), 
short-dashed (green), long-dashed (magenta), dot-dashed 
(red), and dot-long-dashed (yellow) lines present the re- 
sults with N im = 10 49 , 10 50 , 10 51 , 10 52 , and 10 53 s _1 , 
respectively. A Galactic dust-to-gas ratio + filling factor 
(k = 1) is assumed, (b) Right panel: same as Panel (a) 
but with k = 0.1. The regression slopes corresponding to 
D7 1 and n c oc Z> are also presented in the lower- 
left corner in each panel. 



4.4. Results 

In Fig.[3j the results of the static models are plotted over the ob- 
servational samples for various N lon . Hirashita et al. (2001) sug- 
gest that N lon = 3 x 10 49 s _1 on average for Galactic H n regions. 
Indeed N[ on ~ 10 49 -10 5() s _1 is consistent with the Galactic Hn 
region sample as shown in Fig. [3^. In Fig. [3^, the dust-to-gas ra- 
tio is assumed to be Galactic (k = 1), while in Fig.[3]5, k = 0.1 to 
take into account the relatively low-metallicity of the HST sam- 
ple (Table [2]). As shown in Fig. [3^, the size-density relation is 
relatively insensitive to the change of A^n- The reason for this is 
described in the last paragraph in Sect. 14.21 since Tsd » 1, the 
increase of rs with Ni 0n is compensated by the decrease of y\, 
and as a result, r; increases only slightly. 

Because of this weak dependence of r; on M on , it is ex- 
tremely difficult to explain the data of some BCDs, unless we 
assume an extremely large N lon . However, the size-density rela- 
tion of the extragalactic sample, especially that of the BCDs, is 
readily explained if we assume a lower dust-to-gas ratio typical 
of the BCD sample (k = 0.1), as shown in Fig. [3}?. For this value 
of dust-to-gas ratio, Tsd is typically S5 1, and rj increases almost 
in proportion to M on ■ 

Because the size-density relation implies a constant ionized 
gas column density, as outlined in the Introduction, we might in- 
stead expect that the data would be consistent with constant Tdi, 
the dust optical depth within the ionization radius. To test this, 
because the dust content is expected to decrease with decreas- 
ing metallicity, we should correct the inferred ionized gas col- 
umn densities for the different metal abundances of the samples. 
Fig.|4]shows the ionized gas column densities multiplied by their 
oxygen abundance relative to solar plotted against region diame- 
ter. This correction assumes that k oc Z/Z , where Z is metallicity 
or oxygen abundance. The horizontal linefl correspond to con- 
stant Tdi of 0.2, 1, 2, 5, and 10. For these calculations of Tdi, 
M on ranged from 10 48 to 10 53 s~', and electron densities (n e ) 
from 10~ 2 to 10 6 ctrT 3 . A large optical depth of Tdi ~ 10 is only 
possible for large ionization radii achieved with high values of 
Mon(~ 10 52 -10 53 s-') and high densities «n e > ~ 10 4 -10 6 cirT 3 ). 
This is why the Tdi = 10 horizontal line is shorter than the oth- 
ers. Conversely, the low optical depth Tdi = 0.2 value can only be 
achieved for low values of N mn and low densities, (M on ~ 10 48 - 
10 50 s" 1 , <n e > ~ 10~ 2 -1 cm" 3 . This is why the T di = 0.2 line lies 
at large diameters. 



In our static model, the corrected column densities correspond to 



4.9 10 2U T di cm- 
Appendix [A}. 



, assuming a Galactic dust-to-gas ratio of 6 x 10 (see 
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Fig. 4. Ionized gas column densities (cm -2 ) vs. diameter (pc). 
The ionized gas column densities have been multiplied by their 
oxygen abundance relative to solar, assuming the Anders & 
Grevesse ( 1989) calibration. Constant Tdi are plotted as horizon- 
tal lines: Tdi = 0.2 (dot-dashed), 1 (solid), 2 (dotted), 5 (short- 
dashed), and 10 (long-dashed). Symbols are as in Fig. [2] 

It is clear from Fig. [4] that the corrected data are also incon- 
sistent with constant Tdi- Moreover, a linear correction of dust 
column density (or optical depth Tdi) for metallicity does not 
align all the samples. Hence, such a correction is apparently ei- 
ther insufficient or wrong, when considering the dust content of 
low-metallicity SF regions, at least in the context of these simple 
models. 

The above implies that the size-density relation is not a se- 
quence of either A^; on or t^. Rather it is probable that the rela- 
tion should be considered with variable Ni on . This means that we 
should assess how Ni m evolves if we insist that the entire sample 
should be explained by a single "unified" sequence. Thus, in the 
following we include the time variation of A^i on into our model, 
and calculate the consequent variations in density and size. 



5. Evolutionary models 

The treatment used for the evolution of an H n region is based on 
our previous paper, Hirashita & Hunt (2006 hereafter HH06), 
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where we treat the evolution of the number of ionizing photons 
emitted per unit time (Mem) under a given star formation history. 
We extend our models to include the effect of grains according to 
A04. Indeed, as shown below, dust absorption of ionizing pho- 
tons significantly reduces the size of H n regions especially for 
compact Hn regions (Sect. ITTl i. We describe our models in the 
following. 

5.1. Star formation rate 

We assume a spherically symmetric uniform SF region with ini- 
tial hydrogen number density, «ho, and available gas mass for the 
star formation, M gas . Following HH06, we relate the star forma- 
tion rate (SFR) with the free-fall time of gas. The free-fall time, 
tg, is evaluated as 



4.35 x 10 7 

t« = — , — yr ■ 

V"H0 



(10) 



Then, the SFR, tf/(t), is estimated with the gas mass divided 
by the free-fall time as 

, 0.230 m ( -^-) ( _J!™ V' 2 fit) M yr 1 , (1 1) 
\O.l/\lO 7 M /UOOcm- 3 / Jy > oy 

where 6sf is the star formation efficiency assumed to be 0.1 
throughout this paper (HH06), M gas is the total gas mass avail- 
able in the SF region, and the dimensionless function f(t) speci- 
fies the functional form of the SFR. We set the zero point of time 
t at the onset of star formation in the SF region. We investigate 
the following functional form for /(f): 



f exp(-f/f g ) if t > 0, 
\ if t < 0, 



(12) 



where we define the gas consumption timescale t s as 

f g = fff/esF- (13) 

We will refer to this form of the SFR as "exponentially decay- 
ing". In order to examine the effects of continuous (not decaying) 
input of ionizing photons, we also adopt another functional form 
for the SFR called "constant SFR": 



specified in Sect. 15.11 We use the fitting formulae of Schaerer 
(2002) for t,„ and Q(m). We adopt Z = (zero metallicity) for 
the stellar properties following HH06. If we adopt solar metal- 
licity, A^i on decreases by a factor of ~ 2 for a given stellar mass, 
but the decrease of Ni on for metal-rich stars can be compensated 
if we assume ~ 2 times larger SFR. In fact, after taking into 
account the different SFR and initial gas mass, our results are 
comparable to those in Starburst99 (Leitherer et al. 1999). 

5.3. Radius of the ionized region 

In order to treat the range of evolution variations from deeply 
embedded H n regions to normal H n regions, it is crucial to in- 
clude pressure-driven expansion of H n regions. We adopt a sim- 
ple analytical approximation based on HH06. Here, we newly 
include the effect of dust, since A04 have shown that the dust ex- 
tinction significantly reduces the radius of compact H n region. 

We divide the growth of an H n region into two stages: the 
first stage is the growth of ionizing front due to the increase of 
ionizing photons, and the second is the pressure-driven expan- 
sion of ionized gas. The expansion speed of ionizing front in 
the first stage is simply estimated by the increasing rate of the 
ionization radius (Eq. |5J. We denote the Stromgren radius un- 
der the initial density as rso, which is estimated by using Eq. ([1]) 
with «h = "ho- Because of dust extinction, the ionization radius 
r\ is reduced by a factor of y;, which is determined by solving 
Eq. ©, where rsa is evaluated by Eq. (0 with a given D (or k). 
The ionization radius r\ and are estimated by Eqs. © and ©, 
respectively. 

Initially, the increase of r\ is caused by the accumulation of 
ionizing stars. Roughly speaking, as long as r\ (the increase rate 
of the ionization radius) is larger than the sound speed of ionized 
gas, C\\ (we assume Cn = 10 km s _1 in this paper), the ionizing 
front propagates before the system responds hydrodynamically. 
Therefore, we neglect the hydrodynamical expansion if t\ > Cn, 
and adopt the fixed density «eo- 

Once < Cn is satisfied, pressure-driven expansion is 
treated. Since the density evolves, we calculate the Stromgren 
radius rs under the current density «h by using Eq. {T). In this 
situation, the growth of the ionizing region is governed by the 
pressure of ionized gas and the luminosity change of the central 
stars has only a minor effect. Therefore, the following equation 
derived by assuming a constant luminosity (A04) approximately 
holds: 



fit) 



1 if?>0, 

if t < or t > U 



We stop the calculation at t — f g in the constant SFR, since after 
this time the effect of gas consumption should modify the density 
structure around the ionizing source. 

5.2. Evolution of ionizing photon luminosity 

The evolution of the number of ionizing photons emitted per unit 
time, Ni on , is calculated by (HH06) 

Mon00= dm dt' Q(m)4>(m)ip(t-t'), (15) 
Jo Jo 

where r,„ is the lifetime of a star with mass m, Q{m) is the num- 
ber of ionizing photons emitted by a star with mass m per unit 
time, and <p(m) is the initial mass function (IMF). We assume a 
Salpeter IMF <p{m) cc m~ 235 (Salpeter |1955l ) with stellar mass 
range of 0.1-100 M Q . The functional form of the SFR if/(t) is 



(c u 


v 2 

| - &. 




' PI 



(14) A =PR 



1 P n 
i i 

pi 



■Tdi) 



-2 



(16) 



where pi and pn are the gas densities outside and inside of the 
ionized region, respectively, and the function g(T<n) is defined as 



S(Tdt) = — (2e~ Tji - 2 + 2r dl - t\ 



4)- 



The density ratio pn/pi is estimated by (A04) 
2 + e T »(2 - 2t„ + t 2 q ) 



Pn 
Pi 



,. ,3/2/ x 3/2 

no \ 



To 



-2 + e T «(2-2Tdi + rL) 



(17) 



(18) 



where no is evaluated by the Stromgren radius under the current 
ionizing photon luminosity and the initial gas density «ho, and 
to is the dust optical depth through no under density p\. Since «h 
is defined in the ionized region, we relate it with pn as 



pn = P'«h«h ■ 



(19) 
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where yume is the gas mass per hydrogen atom. The density in 
the neutral region is assumed to be constant: 

pi = yum H n H o • (20) 

Eqs. ([Tol l and ([T8l are numerically integrated to obtain r\ as a 
function of t . 

The above pressure-driven expansion is treated as long as 
d(y;rs)/df > 0. When the SFR declines significantly, rs begins 
to decrease. Thus, d(y;rs)/df < may be satisfied at a certain 
time. When d(y;rs)/df < 0, = y\r$ is adopted with fixed tin', 
that is, we finish treating the dynamical expansion. 

6. Initial conditions and results 

6.1. Active and passive star formation 

We have argued in previous papers that star formation at low 
metallicities can proceed in two ways: one "active" mode in 
which stars form in dense, compact complexes, and the other 
"passive" mode where star formation occurs over more diffuse 
and extended regions. The size-density relation presented here 
for extragalactic Hn regions naturally lends itself to the ac- 
tive/passive dichotomy (see also Hunt et al. I2003l l. Dense re- 
gions tend also to be compact, while less dense ones are more 
extended. Here we examine these "active" and "passive" cases. 
In HH06, SBS 0335-052 was used as a prototype of the "ac- 
tive" mode, while IZw 18 the "passive" one. This representation 
was based mainly on the radio continuum results for both galax- 
ies (Hunt et al. 120041 Hunt et al. 120051 1. Their linear emission 
measures differ by 3 orders of magnitude, and the resulting den- 
sities by a factor of 10. However, even in the optical, the electron 
densities inferred from the [S n] line ratios differ by a factor of 
5 or so (SBS 0335-052 has n e ~ 400-600cirT 3 , and IZw 18 
n e ~ 100-120cirT 3 , see Tabled The <n e > of SBS 0335-052 
calculated in this paper (Tabled is much lower than that adopted 
in HH06 (n e ~ 400-600 cm 4 ; Izoto v et al. |1999). As m entioned 
in Sect.[2](see also Kennicutt 1984; Rozas et al. 1998), because 
of dumpiness the densities from the [S n] lines tend to be signif- 
icantly higher than the (n e ). The "active" nature of some BCDs 
therefore may not emerge in the optical, at least with the rms 
densities we have adopted here. However, if a BCD can be clas- 
sified as "active", when observed in the radio, it would have a 
compact, dense (or ultra-dense) H n region (see also Sect. 16.21 ). 

Here, for consistency, we examine "active" and "passive" 
cases by adopting similar initial conditions to those in HH06: 
the active case is modelled by the so-called "compact model", 
and the passive one by the "diffuse model". Appropriate initial 
densities are «ho ~ 10 3 — 10 4 cm 4 for the compact model and 
«ho ~ 10 2 cm 4 for the diffuse model (HH06). In addition, we 
also adopt a denser model as a "super-active" class to explain 
the extragalactic radio sample, and call this model the "dense 
model". For all cases, the gas mass is fixed as M gas = 10 7 M , 
which is similar to the value adopted in HH06. We calculate the 
evolution of the size-density relation of ionized region for the 
continuous/burst star-formation histories (Sect. lBTTT i. 

We also examine the dependence on k, by setting k — 0, 
0.1, and 1 for dust-free, dust-poor, and dust-rich (Galactic dust- 
to-gas ratio) cases, respectively. These cases also correspond to 
varying degrees of dust inhomogeneity or gas filling factor, as 
described in Sect. 14.31 The sizes and densities of H n regions with 
a dust content corresponding to k — 0.01 are not significantly 
affected by dust (i.e., the results are similar to the case with k = 
0), except for the dense model as we discuss later in Sect. 17.11 



On the other hand, the effects of dust extinction can be quite 
pronounced for k > 0. 1 . 

6.1 .1 . Dense model 

We adopt n m = 10 5 cirr 3 (f ff = 1 .4x 10 5 yr) and M gas = 10 7 M . 
With these assumptions, the SFR if/(t = 0) = 7.3M yr~ 1 . First 
we show the evolution of N[ on in Fig.|5^ for the exponentially de- 
caying and constant SFRs. In both SFR scenarios, A^; on increases 
as time passes and stars form. 

We present the basic evolutionary behavior of the size and 
density of H n regions, focusing on the effect of dust, which is 
newly incorporated in this paper. In Figs. [6^ and b, we show the 
time variation of the ionization radius and of the density, respec- 
tively, for various k in the exponentially decaying star formation 
history. The cases with k — (no dust) are the same as those in- 
vestigated in HH06. The general behavior of pressure-driven ex- 
pansion is qualitatively similar in all the dust-to-gas ratios, but 
quantitatively different. The ionization radius decreases as the 
dust-to-gas ratio increases, because of the absorption of ionizing 
photons. The density decreases more rapidly in the dustier case, 
since the expansion speed normalized to the ionization radius is 
larger (for a smaller regions, the expansion at a given speed de- 
creases the density more rapidly). The evolution of «h and r\ is 
quantitatively similar. 

6.1.2. Compact model 

We adopt n m = 3 x 10 3 cm 4 (t B = 7.9 x 10 5 yr) and M gas = 
10 7 M , as before. With these values for the compact model, 
the SFR if/(t = 0) = 1.3M yr . First we show the evolution 
of Ni 0B in Fig. \5]p for the exponentially decaying and constant 
SFRs: Up to t ~ 1 Myr, Ni on increases as stars form. For the 
exponentially decaying SFR, it begins to decrease at t ~ 5 Myr 
because of the death of massive stars, while for the constant SFR, 
it becomes asymptotically constant because the death of massive 
stars is compensated by their continuous birth. 

We present the evolutionary behavior of the size and density 
of Hn regions in Figs. [5J; and d for the exponentially decaying 
star formation history. These models show similar behavior to 
the dense models. The increase of r r and the decrease of «h stop 
at a later stage of the evolution because of the exponential decay 
of the SFR. 

The temporal behavior in the constant star formation is very 
similar to that in the exponentially decaying one. The difference 
is that «h and r\ continue to decrease and increase, respectively, 
because the SFR does not decay. 

6.1.3. Diffuse model 

For the initial density and the radius, we assume «ho = 100 cm 4 
(% = 4.4 x 10 6 yr) and M gas = 1O 7 M , as above. In this case, 
for the diffuse model, the SFR if/(t = 0) = 0.2 M yr 1 . The evo- 
lution of N[ on is shown in Fig. [5}; for the exponentially decaying 
and constant SFRs: around 12 Myr, iVj on is, respectively, 6 x 10 52 
s and 8 x 10 52 s _1 . The behavior of N lon in the diffuse models as 
a function of time is qualitatively the same as that in the compact 
models (Sect. 16. lT2t . but with a different SFR. We present the ba- 
sic evolutionary behavior of the size and density of H n regions 
for the diffuse models with various dust-to-gas ratios in Figs. [(J 
and d for the exponentially decaying star formation history. The 
cases with k — (no dust) are the same as those investigated 
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Fig. 5. Time evolution of the number of ionizing 
photons emitted per unit time (iVi n). The dense 
(Panel a), compact (Panel b), and diffuse (Panel 
c) models are adopted. The solid and dotted lines 
in each panel correspond to the exponentially de- 
caying and constant star formation rates, respec- 
tively. 
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Fig. 6. Time evolution of the ionization radius (n: 
Panels a, c, and e) and the hydrogen number den- 
sity (n H : Panels b, d, and f). The dense (Panels 
a and b), compact (Panels c and d), and diffuse 
(Panel c and d) models are adopted with various 
dust-to-gas ratio normalized to the Galactic value, 
k. The solid, dotted, and dashed lines correspond 
to k = 0, 0.1, and 1, respectively. The exponen- 
tially decaying star formation rate is adopted. 



in HH06. For the dependence on the dust-to-gas ratio, the same 
qualitative discussion as in Sect. |6.1.2l holds. 

The increase of r r and the decrease of «h stop for the same 
reason as in the compact models. The time evolution in the con- 
stant star formation history is very similar to those in the expo- 
nentially decaying one, but «h and r\ do not stop decreasing and 
increasing, respectively. 

6.2. Evolutionary tracks on the size-density relation 

In Figs. [7^ and b, we show the time evolution of the ionized re- 
gion on the size-density diagram for the exponentially decaying 
SFR and the constant SFR, respectively. For comparison with 
the data, we show the diameter of the ionized region, D\. 

Di = 2n. (21) 

Below we discuss the results of the three models. 

6.2.1. Specific models 

We observe from Fig.|7]that the dense models reproduce the data 
points of the radio sample. This means that the radio sample 
can be understood as an extension of the "active" mode of star 
formation toward higher density. However, if a Galactic dust-to- 
gas ratio (k — 1) is assumed, the predicted sizes are smaller than 
those of the radio sample. This implies that the effects of dust ex- 
tinction are suppressed because of an intrinsically low dust con- 
tent, a relatively high gas filling factor, or a non-homogeneous 
dust distribution, but in a way that is not strictly related to metal- 
licity. 

Our models are based on the same gas mass, 1O 7 M , in 
the dense, compact, and diffuse models. This implies that the 



star formation activity traced by the radio sample cannot be ne- 
glected compared with that traced by the HST (optical) sample. 
As shown later in Sect. 17.11 the fraction of ionizing photons ab- 
sorbed by dust tends to be larger in denser regions. Thus, al- 
though they are potentially important to the total star formation 
activity, compact and dense H n regions such as we have assem- 
bled in the radio sample would tend to be overlooked because of 
dust extinction. 

Fig. [7] shows that the evolutionary tracks of the compact 
models in the size-density diagram explain the upper part of the 
HST sample with ages of ~ 3 Myr. Thus, the size and density 
of the "active" class can be explained with the initial condi- 
tions of the compact model, while as shown below, those of the 
"passive" class can be reproduced with the less dense initial 
conditions of the diffuse model. This is consistent with the pic- 
ture that "active" and "passive" star formation processes origi- 
nate from dense and diffuse SF regions, respectively (Hirashita 
& Hunt l2004l HH06). 

The different classes are also distinguished by different ef- 
fects of dust extinction. This is shown in Fig. [8] which presents 
the evolution of the dust optical depth, Tdi, for the different 
classes of models. It can be seen that Tdi is not a linear func- 
tion of k; if k increases, Tsa increases with a fixed values of «h 
and TYion (Eq.|7j, but yj decreases (Eq.[3j. As a result, the product 
Tdi = JiTsd is subject to these competing effects. In particular, 
if Tdi ^ 1, this nonlinearity becomes significant. This is why Tdi 
changes only by a factor of 8 under a change of k by two or- 
ders of magnitude in the dense model. In other words, if dust 
extinction dominates, the size-density relation is expected to be 
aligned on a strip with 1 ^ Tdi ~ 10, although we should also 
consider the possibility that some diffuse regions have low dust 
optical depth if k % 0.01 is appropriate. 
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Fig. 7. Relation between the rms electron number density (n e ) and 
the ionization diameter D\ for a) the exponentially decaying SFR 
and b) the constant SFR. The dense, compact, and diffuse models 
are shown by solid, dotted, and dashed lines, respectively, and the 
three classes for each model correspond to k = 0.01 (top), 0.1, and 1 
(bottom). The ages along the model curves are denoted by different 
symbols as marked. The data of the radio sample and the HST sample 
are shown with the same symbols as in Fig. [2] 
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Fig. 8. Evolution of r di (the dust optical depth over 
the ionization radius). The solid, dotted, and dashed 
dotted lines present the results with k = 0.01, 0.1, 
and 1, respectively. The models adopted are a) dense 
model, b) compact model, and c) diffuse model. 
The exponentially decaying SFR is adopted in each 
model. 



6.2.2. Summary of the evolutionary tracks 

By examining the above three cases, we have inferred that very 
compact dense H n regions such as the radio sample, compact 
H ii regions in "active" BCDs, and diffuse H n regions in "pas- 
sive" BCDs are different populations which have emerged be- 
cause of different initial densities. Our models also suggest that 
denser regions suffer more from the effects of dust than less 
dense ones; as shown in Fig. [8] "active" and "passive" star for- 
mation modes are characterized by different dust optical depths. 

Contrary to the expectation at the end of Sect. 14.41 the ra- 
dio and HST extragalactic samples cannot be reproduced with a 
single initial condition. Rather we should consider that we are 
seeing an "envelope" of individual evolutionary tracks of H n re- 
gions with diverse initial densities. The concept of "envelope" 
naturally suggests that we selectively observe the end points of 
the evolution: i.e., we tend to sample "active" SF regions with 
young ages of an order of Myr, while we observe "passive" re- 
gions with 10 Myr or older ages. In this sense, the observed size- 
density relation of extragalactic H n regions can be taken as an 
age sequence, although the relation contains different popula- 
tions starting from different initial densities. The extreme youth 
of "super-active" (extremely dense) regions is also supported by 
the typical rising thermal radio spectrum from free-free gas ab- 
sorption (Cannon & Skillman 120041 HH06; HH06 also argue 
that the initial density is important for the rising spectrum). 

As mentioned before (Sect. 13.21 ), some of the galaxies in the 
HST sample are also present in the radio sample, as (ultra-)dense 
radio Hn regions. Hence, it may not only be a question of evo- 
lution, but rather also of what the observations are probing. It 
seems that the more finely some Hn regions are probed, the 
smaller and more dense they appear, even in the same source. 
The H ii region in NGC 5253 has a density of 7500-40000 ctrT 3 
if observed with a radio interferometer, but ~ 300ctrT 3 with 
an optical spectrum. The same is true for He 2-10 (3300- 



10000 cm 3 in the radio and ~ 500 cm 3 in the optical; Vacca & 
Conti [19921 IIZw40 (35000 cm- 3 in the radio and ~ 190 cm 3 
in the optical), and, as mentioned above, SBS 0335-052 (1500- 
7900 cm 3 in the radio and ~ 500-600 cm -3 in the optical). 
Thus, we propose that there is a hierarchy of size and density. 
This picture is consistent with that proposed by Kim & Koo 
(|200Tl i, Efremov & Elmegreen (1998) and Elmegreen (2000), 
and will be further discussed in Sect. 17.41 

7. Discussion 

Here we discuss the implications of our results for evaluating 
star-formation activity from ionizing photons, and examine how 
the radio and HST samples fit into the larger picture of star for- 
mation locally and at high redshift. 

7. 1 . The effects of dust 

Even in the dust-poor case with k < 0.1, the absorption of ion- 
izing photons by dust is significant, especially in dense regions 
with «h ^ 10 3 cm~ 3 (see Fig. O. The effect of dust extinction in 
the HST sample is expected to be smaller than in the radio sam- 
ple, but it is never negligible according to our results if k <: 0.1. 
Thus, dust extinction could constrain the size of H 11 regions in 
BCDs, just as it does in Galactic H 11 regions. 

The most straightforward way to confirm if dust extinction 
significantly alters the size-density relation is to overplot on the 
data the dust optical depth calculated by the models in Sect. [4] 
For this purpose, we fix the expression for T& (Eq. [6]( to a con- 
stant, and then use the appropriate model n e and radius (r{) which 
would give that tj; as an asymptotic value. This is done in Fig. 
|9]for r di = 4, 7, and 10 for k = 1 and r di = 1 for for k = 0.01, 
0.1, and 1. This is a similar diagram to Fig. [4] but with the dif- 
ference that there we corrected the ionized gas column densities 
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assuming a linear dependence on metallicity. As for Fig. |H t& 
was calculated for densities which varied from 10~ 2 to 10 6 cm -3 , 
and Nioa from 10 48 to 10 53 s _1 ; this is the reason that the lines of 
constant Tdi do not extend over all radii and densities. 

With a fixed value of k, large n e and/or large r\ are neces- 
sary to increase the dust optical depth t&. Thus, increasing Tdi 
moves the line of constant dust optical depth upward and to the 
right. However, a small value of k with a fixed Tdi shifts the 
line of constant Tdi in the same direction, because for a small 
value of a:, n e and/or r\ need to increase in order to keep t& con- 
stant. Therefore, the line of constant dust optical depth depends 
strongly on the adopted value of k. 

Several salient points emerge from Fig. [9] and from the com- 
parison of Figs. |9j and [4] First, dust is an important factor in most 
of the H ii regions in our sample, both Galactic and extragalactic 
ones, since Tdi k, 1 is indicated. While in Galactic H n regions, 
this is well established (e.g., Gail & Sedlmavr 119791 Melnick 
[19791 Churchwell et al. 119951 , in low-metallicity BCDs the ef- 
fects of dust have generally been thought to be negligible. 

The second point is that in Fig. |9j the data are fairly well 
approximated by different values of constant (except for the 
extremely compact and dense radio sample). Fig. |8] shows that 
Tdi predicted by our evolutionary models varies by no more than 
a factor of two over the gas consumption lifetime. This relative 
constancy of Tdi predicted by our models may be one of the main 
underlying factors in the size-density relation observed in H n 
regions. 

Thirdly, the correction for metal abundance applied in Fig. |4] 
does not align the samples (see Sect. [3). The only difference be- 
tween Fig.|9jand Fig.|4]is the correction for metallicity assuming 
that k oc metallicity. The implication is that this correction is not 
valid, namely that the dust-to-gas ratio is not linear with metal- 
licity (see also Lisenfeld & Ferrara |T998l Hirashita et al. 2002). 

Finally, the dense (radio) sample can only be approximated 
by small values of k ~ 0.01. Because the mean metal abundance 
of the radio sample is ~ 0.22Z o , this low value of k cannot re- 
flect a low dust-to-gas ratio resulting from a linear variation with 
metallicity. Instead, it must point either to regions evacuated of 
dust through stellar winds, or to an intrinsic lack of dust per- 
haps because of extreme youth. In addition, a relatively high gas 
filling factor or an inhomogeneous dust distribution with a cen- 
tral cavity could contribute to the small k. It is likely that all of 
these alternatives are shaping the properties of these extreme H n 
regions. 

The generally high values of Tdi which are consistent with 
our H n-region samples suggest the possibility that the size of 
many H n regions is limited by dust absorption of ionizing pho- 
tons, rather than the consumption of ionizing photons by neutral 
hydrogen. If dust extinction governs the size of the H h region 
(we call this situation "extinction-limited"), we would naturally 
expect a constant column density of H n regions, which is equiv- 
alent to a constant dust optical depth under a fixed mass absorp- 
tion coefficient of dust grains. Natta & Panagia 1 19761 1 also de- 
rive a dust optical depth of order unity in the Lyman continuum 
from Galactic H n regions. 

To further quantify the effects of dust extinction in our mod- 
els, we examine the fraction of ionizing photons absorbed by hy- 
drogen. This fraction is denoted as / and the fraction of ionizing 
photons absorbed by dust becomes 1 -/. By using y\ determined 
in Eq. ©, / is expressed as (Spitzer fl 9781 1 

f = yl (22) 

In Fig. [TpJ we show the evolution of / for the dense, compact, 
and diffuse models. It is clear that more than half of the ionizing 



photons are absorbed by dust even with k = 0.01 in the dense 
model. With k ^ 0.1, the dust absorption is generally severe. 
Thus, we expect that many of the H n regions in BCDs suffer 
significant dust absorption of ionizing photons. This means that 
we would grossly underestimate the total mass of massive stars 
in compact Hn regions in BCDs even from radio observations. 
The effect is even more pronounced for hydrogen recombina- 
tion lines. Such an effect is indeed observed in SBS 0335-052 
(Reines et al. 2008), and more work is needed to establish this 
for statistically significant samples. 

However, as noted in Sect. 14.21 inhomogeneities in the dust 
distribution affect / (Tdi), which may be larger (smaller) if we 
consider, for example, a central dust cavity (Natta & Panagia 
119761 Inoue 120021). If such a cavity forms efficiently on a 
timescale of ~ 10 s yr in compact Hn regions, / may increase 
and the above underestimate of the mass of massive stars could 
be less severe. For further qualitative discussion, one should treat 
a motion equation of dusty gas, which is left for future work. 

In summary, the above estimates of Tdi and / imply that the 
dust grains play a central role in determining the ionization ra- 
dius. In particular, the small values of / in compact ionized re- 
gions indicate that the effect of dust grains significantly reduce 
the number of ionizing photons and the ionization radius. This 
means that there could be a large fraction of massive stars which 
are not traced with radio observations or with hydrogen recom- 
bination lines. 

7.2. Galactic H n regions 

We have shown that the sequence of extragalactic H n regions 
cannot be understood as a sequence with a constant iVj on but that 
it can be interpreted as an "envelope" of various evolutionary 
tracks starting from diverse initial densities. This may also be 
true for Galactic Hn regions. Some Galactic Hn regions are as- 
sociated with a single massive star, but some may be associated 
with a forming stellar cluster, for which we can apply the evo- 
lution model described in Sect. [5] This may be the reason why 
a similar size-density relation ((n e ) 00 D7 l ) fits the Galactic Hn 
regions (Fig.f2](. 

In any case, in agreement with previous work, our models 
with constant Tdi suggest that dust plays an important role in 
Galactic H n regions. As shown in Fig. [9] the more compact H n 
regions are well fit by Tdi ~ 10, while the less compact ones 
by Tdi ~ 2. Hence, the Galactic H n region sequence could be a 
sequence of constant Tdi. Nevertheless, because high values of 
Tdi are only possible with high densities and large Mem, it could 
be that the most obscured regions are also the most intrinsically 
luminous, as suggested above. 

7.3. "Active" and "passive" BCDs 

In Sect. [6] we have shown that it is not possible to reproduce 
the entire sample of extragalactic H n regions with a single ini- 
tial condition. Indeed, the two "active"/"passive" star-formation 
modes proposed by Hunt et al. (120031 ) and elaborated upon in 
subsequent papers (Hirashita & Hunt 2004, HH06), could natu- 
rally result from different initial densities. The difference in the 
gas density (and the compactness) explains the variation in far- 
infrared luminosity and dust extinction (Takeuchi et al. 2003; 
Hirashita & Hunt 120041 Takeuchi et al. 2005 ) and in the molec- 
ular fraction (Hirashita & Hunt 2004) between the two classes. 
We had also hypothesized that the "active" class could be as- 
sociated with the existence of SSCs. This is supported by the 
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Fig. 9. Size-density relation with a constant dust optical 
depth T di . The data symbols are the same as in Fig. [2] (a) 
Left panel: The long-dashed, dot-short-dashed, and dot- 
long-dashed lines show T di = 4, 7, and 10, respectively, 
with k = 1, as predicted by the asymptotic behavior of 
our models, (b) Right panel: The solid, dotted, and dashed 
lines correspond to Tdi = 1 with k = 0.01, 0.1, and 1, 
respectively. 
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Fig. 10. Evolution of / (the frac- 
tion of ionizing photons absorbed 
by hydrogen). The solid, dotted, and 
dashed dotted lines present the re- 
sults with k = 0.01, 0.1, and 1, re- 
spectively. The models adopted are 
a) dense model, b) compact model, 
and c) diffuse model. The exponen- 
tially decaying SFR is adopted in 
each model. 



sizes and densities of the known SSCs in the Antennae galaxies 
(Gilbert & Graham 2007 ), which are well approximated by the 
compact models for the active mode of star formation. 

In this paper, we have added the radio sample as a dense ex- 
treme. The radio sample could be called "super-active" in the 
sense that their H n regions are denser and more compact than 
even the "active" BCDs. Indeed, our models imply that they con- 
stitute a different population from the HST sample, and that they 
are never detected in the optical because of heavy dust extinc- 
tion, even at low values of k. Since a large number of ionizing 
photons are absorbed by dust grains (Sect. 17.11 ). even the lumi- 
nosity of the radio continuum or radio recombination lines un- 
derestimates the total mass of massive stars. 

The "super-active" mode of star formation seen in the ra- 
dio sample could also be consistent with the idea that dense 
conditions foster the formation of SSCs (Billett et al. 20Q2J). In 
fact, ultra-dense H n regions or radio "supernebulae" have been 
proposed as SSC progenitors, observed in their very young em- 
bedded state (Kobulnicky & Johnson 1999, Johnson et al. 2003, 
Johnson et al. 2009). This would be consistent with our find- 
ing that extreme conditions are necessary to explain the radio 
sample. The young age predicted by our models would also be 
consistent with the age of 51 Myr derived from statistical con- 
siderations by Kobulnicky & Johnson ( I1999I ). 

7.4. Density, time, and size in hierarchical star formation 

Elmegreen and collaborators have proposed that there is a 
timescale-size hierarchy in SF regions (Efremov & Elmegreen 
1998 Elmegreen 2000, and references therein). They found a re- 
lation between the time difference (Af) and size and/or separation 
(S) of clouds and star clusters. We have attempted to assess this 
trend for the HST sample in Fig. QT|by plotting the H/J equiv- 
alent width (EW(H/J)) against the diameter of the ELC as de- 
scribed above. After about 3.1 Myr, EW(H/J) linearly decreases 



with starburst age (Leitherer et al. 11 9991 ), so we might expect a 
trend between age and size to be traced by EW(H/3) and ELC di- 
ameter. The figure shows a weak negative correlation, formally 
> 2cr. The two dashed lines, correspond to different methods for 
determining slopes (Isobe et al. [1990): minimizing the perpen- 
dicular distance to the line gives EW(H/?) oc q -0,23 * 0,1 , and the 
ordinary-least-squared (OLS) bisector, EW(H/i) oc £)-0.55±0.i 
The fitted intercepts at zero size correspond roughly to ~ 1000A, 
roughly the EW(Hj6) "plateau" for ages younger than ~ 3.1 Myr 
at oxygen abundances 0.05 Z . The fitted power-law indices for 
the HST sample are very close to that found by Efremov & 
Elmegreen ( flWBl and Elmegreen d2000| ) of Af oc $0.33-0.5 for 
molecular clouds and star clusters in the LMC. Hence, even in 
extragalactic Hn regions, there is some (albeit weak) indication 
that time and size are related, perhaps because of the imprinting 
of turbulent molecular cloud fragmentation on the star clusters 
responsible for ionizing the Hn gas. 

Indeed, the size-density relation of Hn regions, n e oc 
Dr 1 , may be compatible with the above time-size relation. 
Considering that the timescale of star formation scales with the 
free-fall timescale (see Sect. 15. H and Elmegreen 2000), it may be 

— 1/2 

reasonable to assume a scaling relation Af oc n e 1 . Combining 

1 /2 

these two density scaling relations, we obtain Af oc Z> , as long 
the typical size or separation of SF regions is determined by the 
size of ionizing radius, S ~ D\. We would thus obtain the same 
scaling relation as shown by Efremov & Elmegreen ( 19983 an d 
Elmegreen d2000l ). 

Although our data for H n regions cannot be used to deter- 
mine star formation activity from small scales to large (galactic) 
scales (e.g., Elmegreen 200(3), we would argue that pressure- 
driven expansion in H n regions and ionization of the surround- 
ing medium can play a significant role in determining the typical 
extension or separation of SF complexes. Considering that the 
dust optical depth affects the ionization radius of H n regions, it 
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Fig. 11. H/3 equivalent widths (A) vs. Diameter (pc) of the SF 
region in the HST sample. Points are shown as in Fig. Q] The 
regressions correspond to the best-fit lines with slopes of -0.23 
(minimizing perpendicular distance from the line), and -0.55 
(OLS bisector). 



is also possible that dust could be a main factor for determining 
the typical extension or separation of SF regions. 

7.5. Implications for high-redshift star formation 

The number density of gas in high-redshift primeval galaxies is 
theoretically considered to be ^ 10 3 cirT 3 (Norman & Spaans 
119971 1. Therefore, the SF regions in such galaxies may be ex- 
pected to mimic our compact models. This means that a signifi- 
cant fraction of ionizing photons can be absorbed by dust grains. 
According to Fig. [TO] roughly 1/3 of the ionizing photons are 
absorbed in the compact model even when k = 0.01. A pair- 
instability supernova originating from a Population III star is ex- 
pected to supply ~ 1OM of dust grains (Nozawa et al. 2003; 
Schneider et al. 2004]l; if the baryonic mass of the first object is 
~ 10 5 M o (Tegmark et al. 119971 1, the dust-to-gas ratio becomes 



~ 10 (k ~ 10 ). Thus, k ~ 10 is reached relatively quickly, 
and dust extinction becomes important soon after the death of 
the first Population III stars. Our results also suggest that the ra- 
diative transfer of ionizing photons at high z is strongly affected 
by dust extinction. Thus, the cosmic reionization history should 
be reconsidered by taking into account the effects of dust extinc- 
tion. 

It is also interesting to consider the starburst populations at 
lower (but still high) redshifts (z £ 5). Although resolving in- 
dividual dense H n regions is impossible in high-z galaxies, an 
intense starburst at high z would be expected to produce several 
giant Hn regions or massive ELCs. For example, let us assume 
a star-forming region whose gas mass is three orders of magni- 
tude larger than that assumed for the BCDs (i.e., M gas = 10 10 M G ; 
note that this is not the total gas mass in the entire galaxy but the 
gas mass involved in the current star formation episode). With 
«ho ~ 100 crrT 3 (typical density for the diffuse model) we can 
roughly reproduce a typical SFR ~ 200M o yr~' (Eq. ITTb for 
luminous infrared galaxies (Sanders & Mirabel 1996). If the ex- 
ponentially decaying SFR is adopted, this model after t ~ 30 
Myr can also be used to mimic a typical SFR of more mildly 



star-forming populations such as Lyman break galaxies and 
Lya emitters (e.g., Takeuchi & Ishii |2004l Pirzkal et al. [2007l >. 
Moreover, with «ho ~ 10 4 crrT 3 , the SFR roughly becomes a 
typical value for such extreme starbursts (~ 1000M o yr _1 ) as are 
seen in submillimeter galaxies (e.g., Hughes et al. |1998l ). 

Because the Stromgren radius depends only weakly on N mn 
(Eq. the size of H n regions is only one order of magni- 
tude larger even when the SFR {N loa is proportional to SFR) 
increases by three orders of magnitude. Severe dust extinction 
would lessen the increase of H n region size. Nevertheless, this 
implies that high-z starburst population could produce an im- 
mense dense Hn region whose size (^ 100 pc) is an order of 
magnitude larger than the local (dense) extragalactic H n regions. 

In these massive H n regions, if the H n region size becomes 
larger than lOOpc, pressure-driven expansion ceases to be im- 
portant. The reason is as follows: pressure-driven expansion oc- 
curs on a sound-crossing timescale. Since the sound speed in H n 
regions is ~ 10 km s _1 , the sound-crossing time over a lOOpc 
region is 10 Myr, comparable to the typical lifetime of massive 
stars. Thus, such a giant Hn region cannot be sustained long 
enough for the pressure-driven expansion to modify the ionized- 
gas density. Consequently, the density of such a giant H n region 
reflects the initial density averaged over the current H n region 
size. 



7.6. Maximum-intensity starbursts 

The "dense" ("super-active") and "active" star formation modes 
we have modelled far exceed the empirical global star-formation 
intensity limit of ~45 M Q yr _1 kpc 2 found by Meurer et al. 
(119971 ). However, recent observations of submillimeter galax- 
ies and quasar hosts at high redshift have discovered "hyper- 
starbursts" occurring in small regions of ^ l-3kpc in diameter 
(e.g., Tacconi et al. 2006 Walter et al. 2009). These kiloparsec- 
scale starbursts have SFRs per unit area on the order of 100- 
1000 M e yr~' kpc 2 , comparable to those of our models. Similar 
high starburst surface densities are found in Arp 220 but on spa- 
tial scales of a few 100 pc (Scoville et al. fTWTb . With 10 10 M of 
gas, roughly that in Arp 220, these are roughly the spatial extents 
predicted by our evolutionary models for massive starbursts, as 
described in the previous section. 

If dust is a principle factor in shaping the observable prop- 
erties of H ii regions and ELCs, then we might expect shorter 
wavelengths, in particular the ultraviolet, to be unsuitable for 
sampling SFRs in such objects. This would mean that starburst 
intensity limits both locally and at high redshift would need to be 
reassessed at submillimeter or infrared wavelengths where dust 
reprocessing gives a more accurate picture (see also Gao 2008). 

8. Conclusions 

We have investigated the size-density relation of extragalactic 
Hn regions, focusing on those in BCDs. Motivated by the simi- 
larity of size-density relations of extragalactic H n regions with 
Galactic ones, we have modelled and examined the size-density 
relation of ionized regions by considering the effects of dust, star 
formation history, and pressure-driven expansion of H n regions. 
The results have been compared with several samples spanning 
roughly six orders of magnitude in size and density. We have 
shown that the entire sample set cannot be understood as an 
evolutionary sequence with a single initial condition. Rather, the 
size-density relation reflects a sequence with different initial gas 
densities. Thus, a hierarchical structure of SF regions with vari- 
ous densities is implied. 
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We have also found that the size of extragalactic H n regions 
is "dust-extinction limited", in the sense that the dust absorp- 
tion of ionizing photons is significant. This naturally explains 
the observed size-density relation of H n regions as following a 
constant column density of ionized gas, if the dust-to-gas ratio in 
H ii regions is constant. The dust extinction of ionizing photons 
is particularly severe over the entire lifetime of the compact ra- 
dio sample with typical densities of S 10 3 cirT 3 . This means that 
the compact radio sample constitutes a different population from 
the optical samples and that star formation activity in such dense 
regions would be underestimated or missed entirely if we use the 
emission from H n regions (hydrogen recombination lines, free- 
free continuum) as the sole indicators of star formation rate. 
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Appendix A: Simple derivation of Eq. © 

Eq. (0 is derived in Hirashita et al. (2001) based on Spitzer 
dl978l l. Here, we present a simpler derivation of Eq. (0. 

We start from the optical depth at the threshold wavelength 
of the Lyman continuum (912 A) denoted as t^: 

t l - Aiooo/1.086 = 12.0£ B _v, (A.l) 

where tl is approximated with the extinction at A — 1000 A, 
which is related with the B — V color excess, Eb-v, by assuming 
the Galactic extinction curve (Spitzer 1978 ). The factor 1.086 is 
required for the conversion from Aiooo (the extinction at 1000 A 
in units of magnitude) to tl. Assuming that the dust-to-gas ratio 
is the same between the ionized and neutral media (Sect. |4.2l i, we 
simply use the proportionality between the color excess and the 
hydrogen column density derived for the Galactic environment 
(Spitzer 1978 ); i.e., we adopt 

N H = 5.9 x 10 21 £ B _y mag" 1 cm" 2 , (A.2) 
for solar metallicity. Then we obtain 

t l = 2.03 x 10" 21 /Vh ■ (A.3) 

Next, we estimate rsd (the dust optical depth for the ionizing 
photons over a path length equal to the Stromgren radius) intro- 
duced in Eq. ©. Noting that the hydrogen column density over 
the Stromgren radius is given by Nh = "h^s, we obtain tl = Tsd 
as shown in Eq. (0 after using the expression in Eq. (0. If we 
consider the optical depth over the ionization radius, tl = tj; 
should be applied. 
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Table 1. Radio sample 



Name 


d 


n e 


D 


12 + log(0/H) 


ref> b 




(Mpc) 


(cm- 3 ) 


(pc) 






He 2-10 


10.4 


10000 


6.4 


8.06 


1,9 






3300 


6.4 




1 


IZw 18 


14.6 


4 


626 


7.17 


3, 10 


IIZw40 


11.1 


35000 


1.1 


8.13 


3, 10 


M33 


0.84 c 


3700 


3.7 


8.78 


5, 17 


Mrk8 


52.5 


20000 


3.7 


8.51 


4, 13 


Mrk33 


25.5 


9000 


1.9 


8.40 


4, 14 


Mrk 1089 


54.6 


20000 


3.9 


8.04 


4, 12 


Mrk 1236 


29.6 


20000 


0.5 


8.07 


4, 16 


NGC4214n 


3.8 


70000 


0.4 


8.22 


4, 11 


NGC5253 


2.5 


7500 


2.7 


8.17 


1,8 






40000 


0.9 




2 


NGC 6946 


5.6 


3700 


9.9 


9.12 


5, 15 


NGC253 


3.4 


10000 


3.3 


8.88 


5, 17 


Pox 4 


52.5 


7000 


2.6 


7.95 


4, 14 


SBS 0335-052 


53.7 


1500 


22.7 


7.30 


6, 10 


SBS 0335-052-1 


53.7 


4300 


9.6 


7.30 


7, 10 


SBS 0335-052-2 


53.7 


7900 


5.4 


7.30 


7, 10 


Tol 35 


25.2 


20000 


0.5 


8.11 


4, 15 


VIIZwl9 


67.4 


5000 


4.8 


8.71 


4, 14 



a References for the radio interferometric observations: 1) Mohan et al. ( 120011 ); 2) Turner et al. (2000); 3) Beck et al. {2002); 4) Beck et al. 

12DTO ; 5) Johnson et al. 120011 ; 6) Hunt et al. (2004); 7) Johnso n et al. ( [20091 . 
b Reference s for t he oxygen abunda nce: 8) Walsh & Roy |1989); 9) Vacca & Conti JT992T); 10) Thuan & Izotov 120051; 11) Kobulnicky & 

Ski llman {T9961 ; 12) G useva et al. ( [20001 ; 13) Shi et al. ( [20051 ; 14) Kunth & Joubert ( 1985 1; 15) Kobulnicky et al. ( fT999l ; 16) Buckalew et 

al. 120051 : 17) Garnett 120021 . 
c The distance of M33 is taken from Freedman et al. U99U . 
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Table 2. HST sample 



Name 


Dist. " 


12+ 




«e([Sll])° 


<«c> L ' 


Fnp 1 


EW(H/3) s 


Filter 


D" 


D' 


Ref. 1 




(Mpc) 


log(0/H) b 


(X10 4 K) 


(cm- 3 ) 


(cm- 3 ) 


(xlO- 14 
ergcrrT 2 s -1 ) 


(A) 




(pixel) 


(pc) 




to 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


Haro3 


18.5 


8.34 


1.019 


180 


37 


26.4 


146 


F606W 


1 1.6 


70.2 


1,2 


He 2- 10- A 


10.4 


8.34 


1.120 


383 


450 


261.0 


34 


F814W 


12.1 


91.4 


3,4,5 


HS 0822+3542 


12.6 


7.46 


2.000 


390 


7 


0.5 


272 


F814W 


30.2 


69.1 


2, 6, 7, 8 


IIZw40 


11.1 


8.13 


1.302 


190 


77 


8.5 


259 


F814W 


18.0 


36.4 


2,9 


IZwl8-NW 


14.6 


7.17 


1.997 


120 


10 


3.2 


81 


F814W 


35.1 


167.8 


2, 10 


Mrk59-1 


13.4 


8.01 


1.354 


79 


13 


15.1 


150 


F606W 


27.6 


268.7 


2, 11 


Mrk71-A 


4.4 


7.89 


1.568 


150 


163 


87.3 


357 


F814W 


7.5 


23.9 


10, 12 


Mrk209 


5.2 


7.81 


1.561 


100 


57 


14.6 


224 


F814W 


13.4 


25.5 


2, 10 


Mrk930 


74.6 


8.06 


1.230 


58 


6 


6.9 


93 


F606W 


16.7 


905.9 


13 


Mrk996 


21.7 


8.00 


1.500 


450 


108 


20.2 


129 


F791W 


6.0 


95.2 


14 


NGC1140 


19.8 


8.28 


1.024 


70 


6 


6.7 


51 


F814W 


24.5 


352.9 


1, 9, 15 


NGC1156 


6.5 


8.39 


1.630 


30 


353 


240.0 


87 


F814W 


10.7 


12.7 


12, 16 


NGC 1741 -A 


54.6 


8.12 


1.130 


120 


88 


88.9 


137 


F814W 


8.5 


337.5 


3, 9, 15 


NGC 5253-1 


2.5 


8.17 


1.190 


290 


1603 


223.9 


181 


F814W 


8.3 


7.6 


4, 17 


Pox 186 


17.5 


7.74 


1.694 


350 


36 


3.5 


375 


F814W 


4.2 


53.8 


18 


SBS 0335-052 


53.7 


7.30 


2.000 


412 


10 


1.9 


230 


F791W 


3.5 


134.4 


2, 15, 19, 20 


SBS 1415+437 


12.8 


7.61 


1.700 


60 


22 


4.5 


166 


F791W 


8.0 


74.7 


22 


Tol 1214-277 


115.6 


7.55 


1.979 


325 


7 


1.9 


290 


F814W 


4.3 


358.2 


15, 23 


Tol 1924-4 16a 


42.4 


7.94 


1.390 


78 


64 


9.7 


59 


F814W 


2.0 


61.3 


24, 25, 26 


Tol 65 


32.0 


7.54 


1.732 


120 


17 


2.6 


221 


F814W 


4.3 


98.5 


23, 27 


UGC4483 


5.1 


7.58 


1.700 


100 


21 


2.2 


169 


F814W 


12.4 


46.4 


15,21 


UM311 


22.6 


8.31 


0.970 


44 


26 


10.5 


279 


F814W 


2.7 


43.9 


13, 15 


VII Zw 403 


1.5 


7.73 


1.540 


41 


45 


4.7 


153 


F814W 


12.2 


13.1 


10, 12, 15 



a Distance from NED. 

b Oxygen abundance in ionized gas. 

c Electron temperature. 

d Electron number density « e ([S il]) measured from the ratio of AA67 17,6731 A. 

e RMS electron number density (;i c ) inferred from the H/? emission measure with diameters as in Col. 11, corrected for extinction and ionized 

helium as described in the text. 
f H/3 flux in the spectroscopic slit. 
g Equivalent width of H/? emission. 

h Diameters of brightest star-forming complex as described in the text, without the enlargement factor. 

1 Diameters of brightest star-forming complex (1.5x that measured in the continuum) as described in the text. 

J Spectroscopic refere nces: 1) Izotov & Thuan {2004b); 2) Thu an & I zotov ([2005}; 3) Vacca & Conti ( f!992) ; 4) Ko bulnicky et al. ( [19991; 5 ) 
Johnson et al. ( [20001; 6) Kniazev et al. f2000b ; 7) Corbin et al. ( 12005) ; 8) Izotov et al. { 2006 ); 9) Guseva et al. jjOQO) ; 10) Izotov et al. (flWTl ; 
11) Noeske et al. ( |20T)0); 12) H unter & Hoffman (H999) ; 13) Izotov & Th uan ([1998) ; 14) Thuan et al. (fl996) ; 15) Izotov et al. ( 12007) : 16) 
Moustakas & Kennicutt ( |2006T l ; 17) W alsh & Roy ( [1989); 18) G useva et al. ( 12004); 19) I zotov et al. ( fl999) ; 20) Izoto v et al. ( |2006l ; 21) Izotov 
et al. ( [1994) ; 22) Guseva et al. ( 12053) ; 23) Iz otov e t al. J2001bt : 24) Guseva et al. ( 12057) ; 25) Papaderos et al. ( 12006) ; 26) Kehrig et al. ( 12006) ; 
27) Izotov et al. d2004ak 28) Kniazev et al. J200T) . 



